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Abstract 

Superconductivity in the recently proposed ground-state structures of atomic metallic hydrogen 
is investigated over the pressure range 500 GPa to 3.5 TPa. Near molecular dissociation, the 
electron-phonon coupling A and renormalized Coulomb repulsion are similar to the molecular 
phase. A continuous increase in the critical temperature T c with pressure is therefore expected, to 
~356K near 500 GPa. As the atomic phase stabilizes with increasing pressure, A increases, causing 
T c to approach ~481K near 700 GPa. At the first atomic-atomic structural phase transformation 
(~1 - 1.5 TPa), a discontinuous jump in A occurs, causing a significant increase in T c of up to 
764K. 

PACS numbers: 74.20.Pq, 74.10,+v, 74.62.Fj, 74.20.Fg 
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I. INTRODUCTION 



At relatively low pressure, hydrogen exists in an insulating molecular phase. In 1935, 
Wigner and Huntington predicted that sufficient pressure would cause both a molecular-to- 
atomic transition and metallization^. Recent ab initio calculations support these predictions, 
and have revealed the precise details associated with both effects. Calculations based on ab 
initio random structure searching by Pickard and Needs^ as well as McMahon and Ceperley^ 
suggest that the molecular-to-atomic transition occurs near 500 GPa, the latter study also 
revealing a profusion of structures that atomic hydrogen adopts; and exact-exchange calcu- 
lations based on density-functional theory (DFT) by Stadele and MartirP suggest a metal- 
lization pressure of at least 400 GPa. In 1968, Ashcroft predicted an even further transition 
in high-pressure hydrogen, a metallic-to-superconducting one~ Within the framework of 
Bardeen-Cooper-Schrieffer (BCS) theory^, three key arguments support this prediction: (i) 
the ions in the system are single protons, and their small masses cause the vibrational en- 
ergy scale of the phonons to be remarkably high (e.g., k^(u) ~ 2300K near 500 GPa - see 
below), as is thus the prefactor in the expression for the critical temperature T c ; (ii) since 
the electron-ion interaction is due to the bare Coulomb attraction, the electron-phonon 
coupling should be strong; and (iii) at the high pressures at and above metallization, the 
electronic density of states iV(0) at the Fermi surface should be large and the Coulomb 
repulsion between electrons should be low, typical features of high-density systems. These 
arguments will be revisited, and demonstrated to indeed be the case, below. 

Ever since the prediction of high-T c superconductivity in hydrogerP, a large number of 
efforts have focused on determining the precise value(s) of tJ^S. In the molecular phase, 
the high-pressure metallic Cmca structure (which transitions to the atomic phase 2 ' 5 ) has 
recently been studied in-depttP^^, and shown to have a T c that increases up to 242K near 
450 GPa. In the atomic phase, estimations of T c have varied widely, but in general suggest 
a large increase with pressure^^. Early estimates suggested that T c « 135 - 170K near 400 
GPa (although, it is now believed that this is inside the molecular phase^, as discussed 
above)^, near 480 - 802 GPa, more recent estimations suggest that T c m 282 - 291K^ : ; and 
near 2 TPa, calculations suggest that T c can reach ~600 - 631K in the face-centered cubic 
(fee) lattice^^. The latter two studies will be discussed further below. 

However, previous studies of superconductivity in the atomic phase have simply assumed 
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FIG. 1: (color online). Ground-state structures of atomic metallic hydrogen, (left) Conventional 
unit-cell of IA\/amd at 700 GPa. (right) 2x2x1 supercell of R-3m at 2 TPa. a and c parameters 
are shown in the figure, as discussed in the text. Fictitious bonds have been drawn for clarity. 

candidate ground-state structures, in a number of cases the fee lattice™™™ 3 . Recently 
though, McMahon and Ceperley demonstrated that such structures are incorrect, and pro- 
vided a comprehensive picture of the (presumably correct) ground-state structures from 500 
GPa to 5 TPeP. Molecular hydrogen was shown to dissociate near 500 GPa, consistent with 
the predictions of Pickard and Needs^. With increasing pressure, atomic hydrogen passes 
through two ground-state structural phases before transforming to a close-packed lattice, 
such as fee or possibly the hexagonal close-packed (hep) lattice. The first is a body-centered 
tetragonal structure with space-group IA\/amd (Hermann-Mauguin space-group symbol, 
international notation) with a c/a ratio greater than unity, as shown in Fig. [I] Including es- 
timates of proton zero-point energies (ZPEs), IAi/amd was demonstrated to transform into 
a layered structure with space-group R-3m near 1 TPa, also shown in Fig. [TJ which is similar 
to a possible high-pressure phase of lithiuirPH R-3m remains stable to ~3.5 TPa, compress- 
ing to a close packed lattice. Given such novel crystal phases and that T c can be very 
sensitive to structural details^, and that modern methods of calculating values of T c should 
be more accurate than those used in earlier studies, it is of great interest to re-investigate 
the long-outstanding predictions of superconductivity in atomic metallic hydrogen. 

This Article is outlined as follows. In Section [TTJ the theoretical background used for 



estimating T c (in this work) is presented. Computational details are given in Section III 
In Section [TV] properties of the ground-state structures of atomic metallic hydrogen as a 
function of pressure, such as lattice parameters and vibrational properties influencing the 
IAi/amd — > R-3m transition, are presented and discussed. Superconductivity is investigated 
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in Section [V] Section fVTI concludes. 



II. THEORETICAL BACKGROUND 



(1) 



According to the BCS theory of superconductivity^, there is a simple relationship between 
T c , the average phonon frequency (a;), N(0), and the pairing potential V arising from the 
electron-phonon interaction, 

feT c = 1.14( W >exp [-j^y 

where k-Q is Boltzmann's constant. This relation is valid as long as /cbT c <C (oo), correspond- 
ing to weak coupling - see below. 

McMillan later solved the finite-temperature Eliashberg equations for Tj^, which includ- 
ing a correction by Dynes^ can be written as 



M 

1.2 



1.04(1 + A) 



GXP [ A-/i*(l + 0.62A)J ' (2) 
where A is the attractive electron-phonon-induced interaction and /i* is the renormalized 
Coulomb repulsion. In high-density atomic hydrogen, Ashcroftf^ demonstrated via an ab 
initio calculation that fi* = 0.089, which is remarkably close to //* = 0.085 obtained from 
the Bennemann-Garland formula^, both results similar to the somewhat standard value for 
high density systems of fi* ~ 0.1. In this work, we therefore take fi* = 0.089 for estimating 
T c . It should be noted that this approximation fails in molecular hydrogen 15 , as investigated 
thoroughly in Refs. [20H221 using a specialized formulation of DFT for superconductivity 
where /i* is calculated ab initio. Interestingly though, at high densities (e.g., near molecular 
dissociation) fi* is found to nonetheless be 0.08 for pressures just above 460 GPeP. 

For A > 1.3 (which in fact corresponds to the situations considered herein), Eq. ^ often 
provides a lower bound to T c . In this case, both a strong-coupling correction as well as a 
correction for the shape-dependence of T c with (u) must be made. These corrections will 
be shown to be especially important in atomic metallic hydrogen, where both A and (oj) are 
large. These corrections are both included in the Allen-Dynes equatio^*' 



k-p,T r 



j, n Win 

i~2 exp 



1.04(1 + A) 



A-/i*(l + 0.62A)_ 
where u\ n is the logarithmic average frequency [i.e., ln(a;i n ) = (In a;) ] and 

1/3 



h 



1 + (A/AO 372 



(3) 



(4) 
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/2 = 1 + X2 , A2 ( 5 ) 



- 1 ) a 2 

A 2 + A 2 

denote the strong-coupling and shape corrections, respectively, where a) 2 = (w 2 ) 1 ^ 2 and Ai 
and A 2 are fitting parameters (e.g., to full solutions of the Eliashberg equations). 
In the original Allen-Dynes equation^, 

Ax = 2.46 (1 + 3.8/2*) (6) 

A 2 = 1.82 (1 + 6.3//) (u> 2 /ou ln ) . (7) 

However, a least-squares analysis between T c as predicted by Eq. ^ and that calculated 
numerically in the Eliashberg formalism for an fee lattice of atomic metallic hydrogen at 2 
TPsP^ suggests the following reparametrization 

Ai = 2.26(1 - 1.28//*) (8) 

A 2 = 2.76(1 + 8.86/2*) (u 2 /u in ) , (9) 

which interestingly provides more accurate values of T c for a selection of low-temperature su- 
perconductors as welW. In passing, we note that there is a very recent further reparametriza- 
tion by the same authors^ that appears especially well-suited for calculating T c for a range 
of /t* values (which could be useful for studying both the molecular and atomic phases 
concurrently, for example). 

In this work, estimates of T c are made using both Eqs. ^ and (J3J) as well as both 
parametrization for Ai and A 2 to give a range of values for T c . 



III. COMPUTATIONAL DETAILS 

All calculations were performed using the Quantum ESPRESSO ab initio DFT cod<P. 
A norm- conserving Troullier-Martins pseudopotential 28 with a core radius of 0.65 a.u. was 
used to replace the 1/r Coulomb potential of hydrogen, which is sufficiently small to ensure 
no core-overlap up to the highest pressure considered in this work (3.5 TPa). The Perdew- 
Burke-Ernzerhof exchange and correlation functional was also used for all calculations, as 
was a basis set of plane- waves with a cutoff of 120 Ry, giving a convergence in energy to better 
than ~0.2 mRy/proton, and 24 3 k-points for Brillouin-zone (BZ) sampling. Phonons were 
calculated using density functional perturbation theory as implemented within Quantum 
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ESPRESSO. Additional computational details pertaining to calculations of the electron- 
phonon interactions will be provided and discussed in Section |V} 

IV. GROUND-STATE STRUCTURES OF ATOMIC METALLIC HYDROGEN 

In this section, we discuss the structural changes that occur in atomic metallic hydrogen 
as a function of pressure. On the basis of our previous study^, we consider I^/amd at 
pressures from 500 GPa to 1.5 TPa and R-3m from 1 to 3.5 TPa. We first consider the 
lattice changes that occur (e.g., compression). We then consider the IAi/amd — > R-3m 
transition and discuss the vibrational properties of each structure that contribute to it, in 
anticipation of the results that are to follow in Section [Vj A further discussion of the ground- 
state and metastable structures of atomic metallic hydrogen can be found in Ref. E] and a 
thorough discussion of tetragonal structures of atomic hydrogen (including the IAi/amd 
structure) can be found in Ref. HSU 

A. Lattice Parameters 

In terms of their primitive unit-cells, I4i/amd is tetragonal with a = b ^ c with two 
symmetry inequivalent atoms at Wyckoff positions (0,0,1/2) and (0,1/2,3/4), and R-3m 
is hexagonal (also with a = b ^ c) and a single symmetry inequivalent atom at the origin. 
The lattice parameters of both structures can therefore be specified completely by a and 
the c/a ratio, as indicated in Fig. [I] For the pressure ranges under consideration, the 
lattice parameters and corresponding Wigner-Seitz radii r s are shown in Tables [I] and [TTJ, 
respectively. 

Between 500 - 700 GPa, IAi/amd resists compression along the c axis, as can be seen in 
the c/a ratio which increases from 2.545 to 2.764. Above 700 GPa the resistance continues, 
but the compression becomes much more uniform. For example, by 1.5 TPa the c/a ratio 
increases to only 2.849. In R-3m, on the other hand, the c/a ratio remains relatively constant 
near 3.05 - 3.06. However, near the predicted transition pressures of ~1 and 3.5 TPa (see 
below and Ref. [3]) there is a preferred compression along the c axis. In fact, not including 
the ZPE suggests that R-3m continues to compress along the c axis to fee above 5 TPa 3 . 
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Pressure (TPa) 


a (a.u.j 


c/a 


r s (a.u.j 


0.5 


2.299 


2.545 


1.226 


0.6 


2.227 


2.599 


1.197 


0.7 


2.134 


2.764 


1.170 


0.8 


2.094 


2.769 


1.149 


0.9 


2.058 


2.774 


1.130 


1.0 


2.027 


2.778 


1.113 



1.5 1.893 2.849 1.049 



TABLE II: Lattice parameters 


and corresponding 


Wigner-Seitz radii r s 


of R-3m as a function of 


pressure. 








Pressure (TPa) 


a (a.u.) 


c/a 


r s (a.u.) 


1.0 


1.832 


3.236 


1.111 


1.5 


1.758 


3.061 


1.047 


2.0 


1.685 


3.054 


1.002 


2.5 


1.629 


3.051 


0.969 


3.0 


1.584 


3.047 


0.942 


3.5 


1.564 


2.943 


0.919 



B. I4i/amd — >■ R-3m Transition 

Static-lattice enthalpy calculations indicate that I4i/amd transforms to R-3m near 2.5 
TPa, but dynamic-lattice calculations (in the harmonic approximation) suggest that this 
pressure is significantly reduced to ~ 1 TPaP. In this section, we use the harmonic and 
quasiharmonic approximations to further investigate the I&i/amd — >■ R-3m transition, in 
anticipation of the results that are to follow in Section |VJ 

Ground-state enthalpies for IA\/amd and R-3m (defined by the parameters in Tables [T] 



and lip were calculated at 1 and 1.5 TPa; Table III ZPEs at each pressure were estimated 
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TABLE III: Ground-state enthalpies and zero-point energies of IAi/amd and R-3m at 1 and 1.5 
TPa. Pressures P are in TPa and enthalpies and energies are in Ry/proton. 





P = 1.0 


1.5 


H 


IA\/amd 


-0.49955 


-0.32022 


R-3m 


-0.49534 


-0.31768 


Ezpe 


IA\/amd 


0.02708 


0.03120 


R-3m 


0.02395 


0.02769 


H + Ezpe 


IA\/amd 


-0.47247 


-0.28902 


R-3m 


-0.47140 


-0.28999 



using the harmonic approximation: Ezpe — J du F(u)hu/2, where F(u) is the phonon 
density of states (PHDOS), and are shown in Table |III| as well. Neglecting zero-point 
pressures and making the simple approximation that the total enthalpies are given by H + 
Ezpe (as was done in Ref. |3j) suggests that the IAi/amd — > R-3m transition occurs nearly 



midway between 1 and 1.5 TPa (see Table III), which is very close to, but slightly higher 



than our original estimate of 1 TPa 3 . Going beyond this approximation, the total enthalpies, 
including the zero-point pressures, can be estimated using a linear approximation, 

-f^tot = H avs + Ezpe, avg + PzPEKvg (10) 

where 

OEzpe 

PZPE= gy- (11) 

is the zero-point pressure, V is the volume, and the subscripts avg denote the average values 



of each quantity between 1 and 1.5 TPa. Note that the latter two quantities in Eq. (10) 
correspond to the zero-point enthalpy. Estimating pzpe using a simple finite-difference gives 
total enthalpies of —0.35765 and —0.35976 Ry/proton for I4i/amd and R-3m, respectively. 
This suggests that the actual transition pressure is a bit lower than the simple enthalpy 
estimate, and is in fact in agreement with our original prediction of ~1 TPa 3 . 
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FIG. 2: (color online). PHDOS of IAi/amd (dashed red line) and R-3m (solid blue line) at 1.5 
TPa. The results have been normalized by the number of atoms per primitive unit-cell. Negative 
values indicate imaginary frequencies. 



As can be inferred from Table |III| and the discussion above, the large decrease in the 
IAi/amd — > R-3m transition pressure from the static-lattice prediction (~2.5 TPa 3 ) arises 
primarily from the significantly lower -Ezpe m R-3m, as well as a more minor contribution 
from the lower Pzpe- To help understand this, the PHDOS for both structures is shown 
in Fig. [2j It can be seen that there are three differences that lead to this behavior: (i) 
the density of high-frequency phonons is greater for IAi/amd, and also occurs at higher 
frequencies (3180 - 4430 cm -1 vs 3000 - 4230 cm -1 ); (ii) IAi/amd has a significant density 
of mid-frequency phonons (~1400 - 3000 cm -1 ), while such modes are mostly absent in 
R-3m (e.g., IA^/amd shows significant peaks at 1510, 2150, and 2990 cm -1 ); and therefore 
(iii) the PHDOS for R-3m is mostly concentrated at low frequencies (< 1400 cm -1 ). 

In passing, we note that R-3m shows a small density of imaginary phonon states at 1.5 
TPa. However, estimating the resulting energy within the harmonic approximatiorP shows 
that it only integrates to 1.372 • 10 -5 Ry/proton. While this is within the accuracy of our 
calculations, this behavior is in fact expected considering that it is indicative of instability in 
a lattice of ions treated classically; and classically, the I4i/amd — > R-3m transition occurs 
near 2.5 TPaP, as discussed above. This is further confirmed by the fact that the instability 
goes to zero with increasing pressure, while such behavior begins to develop in IAi/amd - 
see Ref . [3] 
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FIG. 3: (color online). Gibbs free energy vs temperature for IA\/amd (dashed red line) and R-3m 
(solid blue line) at 1.5 TPa. 

Considering that the PHDOSs are quite different between IAi/amd and R-3m and it is 
finite-temperature effects that are focused on below (i.e., T c ), the possibility of vibrational 
entropic stabilization of one phase over the other exists. In order to estimate this, the 
quasiharmonic approximation can be used, 

dus F(uj) In sinh I - 

'o 

where F(V,T) is the Helmholtz free energy at volume V and temperature T and Eq(V) is 
the static-lattice energy. From this, the Gibbs free energy G can be calculated via G = 
F + pV, given the pressure p. At T = OK, p is given by the external pressure plus the 



POO 

F(V, T) = E (V) + k B T / duo F(u) In 

Jo 



(12) 



zero-point pressure [Eq. (11)]. However, for a fixed V, p is actually a function of T, due to 
thermal expansion of the lattice caused by anharmonic phonons. Contrary to the expectation 
that such effects may be large^, calculations of the melting line of hydrogen (not shown )P^ 
indicate that in fact thermal expansion is in fact small, at least up to a few hundred K where 
atomic metallic hydrogen is likely to melt anyway; and since the purpose of this discussion 
is just to understand qualitative changes that may arise at finite-T, we can estimate p using 
the T = OK value. 



Figure 3 shows the resulting free energy G estimated using Eq. ( 12 ) and the value of 



p at T = OK. Despite the remarkably different PHDOSs (see again Fig. |2J), the behaviors 
of G with T are rather similar for both structures. Thus, temperature is not expected to 
significantly affect the I^/amd — >■ R-3m transition. 
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Based on these results, below we consider IAi/amd from 500 GPa to 1.5 TPa and R-3m 
from 1 to 3.5 TPa, and the IAi/amd — > R-3m transition to occur between 1 and 1.5 TPa. 



V. SUPERCONDUCTIVITY 



In this section, we investigate superconductivity in the IAi/amd and R-3m structures of 
atomic metallic hydrogen. We first provide relevant computational details not discussed in 
Section III, and convergence of the parameters necessary to evaluate Eqs. (|2| and ^ with 



respect to q-points. We then present and discuss the calculated parameters, and finally use 
them to calculate T c values. 

A. Computational Details 

In order to estimate T c using Eqs. ^ and (to), tu\ n , C02, and A must all be deter- 
mined. Of course, the frequency parameters can be calculated directly from the PHDOS. 
For example, (u>) = (l/n^) J duo F(u)u, where n p h is the number of phonon modes and 
J duj F(lu) = n ph . In order to calculate A, a (slowly convergent) double-delta integration 
must be performed on the Fermi surface - see Ref. [33] for a complete discussion and the 
precise implementation details within Quantum ESPRESSO. In order to accurately perform 
this integration, very dense k-point and q-point grids must be used. For the IAi/amd and 
R-3m structures, we found that an electronic grid of 48 3 k-points (and using 24 3 k-points to 



calculate phonons, as discussed in Section III) gave convergence with no discernible error. 



In order to determine a sufficient density for the q-point grid, we performed a series 
of calculations with l 3 , 2 3 , 4 3 , 6 3 , and 8 3 q-points, using IAi/amd at 500 GPa as a test 
case (we also considered R-3m at 2 TPa - not shown). It should be kept in mind that such 
rigorous testing with respect to q-points is especially important in atomic metallic hydrogen, 
as inadequate sampling has been shown to cause significantly incorrect results^ 3 ^. In fact, 
our calculations below suggest that the results of a recent study considering Cs-IV (which 
also has the IAi/amd structure) over a more narrow pressure range than considered here 
using only 3 3 q-points^ gives somewhat incorrect values for A, both in magnitude and trend 
with pressure. The values of A for the various densities of q-points, as well as values of 



lo) and win, are shown in Table IV Relative convergence in A is seen to require at least 
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TABLE IV: Convergence of (u), u)\ n , and A with the number of q-points for IA\/amd at 500 GPa. 



No. of q-points 


(u) (K) 


Win (K) 


A 


l 3 


1660 


1438 


17.91 


2 3 


2307 


1953 


2.82 


4 3 


2277 


2031 


2.06 


6 3 


2287 


1997 


1.67 


8 3 


2295 


2068 


1.81 



6 3 q-points (to be within 10% of the converged value, for example). This is likely due to 
Fermi surface "hot spots" that have been shown to exist in other alkali metals^, which can 



significantly contribute to the electron-phonon interaction. Table IV also shows, on the 
other hand, that (u) and u\ n achieve relative convergence with as little as 2 3 q-points, which 
is consistent with the density found necessary in our previous work to accurately calculate 
the ZPEs of the structures of atomic metallic hydrogerP. Herein, 8 3 q-points were used for 



all calculations (including those in Section IV), corresponding to 59 and 150 total q-points 



in the irreducible BZ for IAi/amd and R-3m, respectively. 



B. Superconducting Parameters 



As shown in Fig. |4j (u) and u\ Q are both extremely high, and increase significantly with 
pressure^ (oo) increases from 2295K to 4056K as the pressure is increased from 500 GPa 
to 3.5 TPa, while coi Q is significantly less (especially for R-3m), increasing from 2068K to 
3308K over the same range. Furthermore, there is a significant decrease in both (a;) and 
wi n at the I^i/amd — > R-3m transition (e.g., by 765K and 926K, respectively, at 1.5 TPa), 
consistent with the results and discussion in Section IIVBI 

More interesting is the behavior of A with pressure; Fig. [5j Near molecular dissociation 
(~500 GPa), the values of A in both the atomic and molecular phases are remarkably close. 



In IAi/amd, A ~ 1.81 (see also Table IV), whereas in the molecular phase (Cmca) A ~ 2 just 
above 460 GPa, but appears to slowly decrease with increasing pressure - see Refs. |2"0"||2"21 
Thus, given that A and fi* are similar in both phases near molecular dissociation (see again 
Section O for a discussion of //*), a smooth variation in T c is likely to occur with increasing 
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FIG. 4: (color online). Temperature prefactors (u) (solid blue line) and uj\ n (dashed red line) as a 
function of pressure in atomic metallic hydrogen. 

3.5 
3.0 
2.5 
2.0 
1.5 
1.0 

0.5 1.0 1.5 2.0 2.5 3.0 3.5 

Pressure (TPa) 

FIG. 5: (color online). Electron-phonon-induced interaction A as a function of pressure in atomic 
metallic hydrogen. 

pressure in this range. 

A large increase in A is seen to occur from 500 - 700 GPa, from 1.81 to 2.32. To help 
understand this, the electron-phonon spectral function, a 2 F(u), at 500 GPa is compared to 
that at 700 GPa in Fig. [6j It can be seen that there is an increase in coupling to both the low- 
and high-frequency phonon modes as the atomic phase stabilizes, while there is relatively 
little change in the coupling to those at mid frequency. The former increase is unexpected, 
as with increasing pressure the PHDOS shifts to higher frequencies, as is indicated in Fig. 
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FIG. 6: (color online). Electron-phonon spectral function a 2 F{uj) of IAi/amd at 500 and 700 GPa. 




-500 



500 



1500 2500 
io (cm 1 ) 



3500 



4500 



FIG. 7: (color online). Electron-phonon spectral function a 2 F(uj) of IAi/amd (dashed red line) 
and R-3m (solid blue line) at 1.5 TPa. 



4j The sharp increase in A, along with the increased (u) and u\ n (see again Fig. |4j), suggests 
that a correspondingly large increase in T c should occur over this small pressure range, which 
is shown below to indeed be the case. 

At the IAi/amd — > R-3m transition near 1.5 TPa, a large jump in A occurs, from 1.43 to 
3.39. This can be understood by comparing a 2 F(u) for both structures; Fig. [Tj In R-3m, 
the large value of A is seen to occur from a strong coupling into the low-frequency modes 
[A = 2 J du a 2 F(u)/u}. This appears to be due to the correspondingly high PHDOS at low 



frequencies, which is absent in IAi/amd (see again Section IV B). Comparing Figs. KHandFf 
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FIG. 8: (color online). Electron-phonon spectral function a 2 F(u) of R-3m at 2 and 3 TPa. 



also shows that in IA\/amd there is decreased coupling into all modes with an increase in 
pressure above 700 GPa, especially at low frequencies. 

With increasing pressure, A in R-3m decreases from its maximum to ~1.98 by 3.5 TPa. 
Figure [8] shows that this results from a weakened coupling into the low-frequency modes that 
was responsible for the sharp increase in A in the first place (near the I4i/amd — > R-3m 
transition). This is likely due to an overall decrease in the PHDOS at low frequencies with 
increasing pressure (not shown). These results, combined with those above, indicates that 
coupling into the low-frequency modes are the key to achieving a large value of A in atomic 
metallic hydrogen. 

The A values presented above are much less than for the (unstable) fee lattice. For 
example, at 2 TPa A ~ 3.06 compared to A ~ 7 - 7.32! 16 * 17 l The large difference can 
be attributed to the even higher PHDOS at low frequencies in fee compared to _R-3?tP, 
which was above suggested to lead to very strong electron-phonon coupling. In Ref. [TBI the 
strong coupling into the low-frequency modes in fee was attributed to the lattice being close 
to instability. While this is consistent with (and likely influential) in the strong coupling 
into R-3m, which is also close to lattice instability near 1-1.5 TPa with protons treated 
classically (see Section IV B and Ref. 3 ), this is not necessarily the cause. For example, such 



behavior does not always occur, as A for IAi/amd appears low near its pressure limits of 
lattice instability (e.g., ~500 GPa), while it becomes largest near the center of this range 
(e.g., ~700 GPa) - see above. 
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FIG. 9: (color online). Values of T c for atomic metallic hydrogen calculated using Eqs. @ and (§]). 
T c for the high-pressure molecular phase is shown using a symbol. 

C. T c 

Using the parameters in Figs. [4]and[5j Eqs. ^ and ^ were used to calculate T c ; Fig. [9j 
The values are seen to be remarkably high, but nonetheless consistent with the discussion 
above. The Allen-Dynes equation and its reparametrization 17: , in most cases, give much 
higher estimates than the McMillan formula (as expected based on the discussion in Section 
[iT] ). Given that uj\ n is significantly less than (co), the increase is thus due entirely to the 
correction factors f\ and f2- Comparing these, in a number of cases, shows that it is fi 
(the strong-coupling correction) that is most important, especially in the reparametrized 
Allen-Dynes equation^. For example, at 700 GPa fx ~ 1.31 and /2 ~ 1.03. 

Just above molecular dissociation, T c ps 315 - 356K. The increase in A combined with 
increases in (u) and u\ n with pressure cause T c to increase up to 403 - 481K by 700K. With 
increasing pressure, T c then decreases (in the IAi/amd phase). However, at the I4i/amd — >■ 
R-3m transition, a large jump in T c then occurs, from 370 - 377K to 561 - 703K. This is 
due entirely to the jump in A, considering that (u) and u\ n are significantly less in R-3m 
(see Fig. |4]) . Although, with increasing pressure, T c then decreases. Thus, ~764K represents 
an approximate upper bound to T c in atomic metallic hydrogen, and possibly conventional 
superconductors (i.e., those described by BCS theory) in general. It is interesting to note 
that secondary maxima in T c occur in both IAi/amd and R-3m. Given that there appears 
to be monotonic decreases in A above their maxima in both structures (see Fig. [5]), this 




••••McMillan 
— — -Allen—Dynes 
^Allen—Dynes (reparam) 



.5 1.0 1.5 2.0 2.5 3.0 3.5 

Pressure (TPa) 
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behavior is simply due to an interplay between this and (u) or u\ n . 

VI. CONCLUSIONS 

In conclusion, we investigated superconductivity in the ground-state structures of atomic 
metallic hydrogen over the range 500 GPa to 3.5 TPa. Near molecular dissociation, the 
electron-phonon coupling A and renormalized Coulomb repulsion in the atomic phase were 
demonstrated to be similar to the values of the molecular phase. This suggests a continuous 
increase in T c with pressure during the molecular-to-atomic transition, to ~356K near 500 
GPa. As the atomic phase stabilizes with increasing pressure, A increases causing T c to 
increase to ~481K near 700 GPa. Near the first atomic-atomic structural phase transfor- 
mation near 1.5 TPa, a large jump in A occurs due to a high PHDOS at low frequencies, 
increasing T c to as high as 764K. 

While the T c values presented incredibly high, they are nonetheless reasonable. However, 
there are two caveats. First of all, even the lowest pressures considered in this work are higher 
than those currently obtainable experimentally (342 GPeP^). Nonetheless, all of them are 
important to planetary physics (albeit likely at temperatures even higher than the values 
of T c ). The other caveat is that it is quite possible that the T c values are higher than the 
melting temperatures of the phases of atomic metallic hydrogen. However, this suggests the 
interesting possibility that the atomic metallic solid phase of hydrogen (at least the I4i/amd 
and R-3m structures) may exist entirely in superconducting states. 
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